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Abstract 

We study the nuclear effects in the electromagnetic structure function F2{x, Q"^) 
in the deep inelastic lepton nucleus scattering process by taking into account 
Fermi motion, binding, pion and rho meson cloud contributions. Calcu- 
lations have been done in a local density approximation using relativistic 
nuclear spectral functions which include nucleon correlations. The ratios 
Rp2{^iQ'^) = 'af'^Ix^q'^) obtained and compared with recent JLab results 
for light nuclei with special attention to the slope of the x distributions. This 
magnitude shows a non trivial A dependence and it is insensitive to possible 
normalization uncertainties. The results have also been compared with some 
of the older experiments using intermediate mass nuclei. 

Keywords: Structure function. Nuclear medium effects. Deep inelastic 
scattering. Local density approximation 



1. Introduction 

Recently Jefferson Lab (JLab) pj] using a high intensity electron beam 
of energy 5.767 GeV has measured the nuclear dependence of the structure 
function in some nuclei by studying the ratio R(x,(5^) = ^^, where a"^ is 
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the inclusive cross section in nuclei and is the inclusive cross section 
in deuterium. The experimental results for the ratio R(x, Q^) have been 
presented by them ll| for 0.3 < x < 0.9 and have re-confirmed the older 
EMC results ji], [s], |J that the structure function of a nucleon is modified 
when it is placed inside a nucleus. This experiment shows that the slope 
of the EMC effect does not scale with the nuclear density and therefore the 
simple models to implement these nuclear effects, based on A or average 
density fits, for example as described by Gomez et al. (sl, fail to describe the 
new and precise results for light nuclei. 

The behavior of Rp2{x,Q^) can be broadly divided into four categories 
viz. X < 0.1 is the shadowing region, 0.1< x <0.3 is the anti-shadowing 
region, 0.3< x <0.8 is the EMC region and beyond x ~ 0.8, known as the 
Fermi motion region. Theoretically, many analysis have been done to study 
the EMC effect and various models have been proposed and discussed in the 
literature jo], 0, Isl, [oj. Several phenomenological parameterizations for the 
nuclear parton distribution functions(NPDFs) have been discussed in the 
literature like the works of Hirai et al. llOl , Eskola et al. Illi , Schienbein et 



al. |12l . |13| which successfully reproduce the nuclear modifications in the deep 
inelastic lepton-nucleus and neutrino-nucleus scattering. 

In this work, we study the nuclear medium effects on the structure func- 



tion within a model based on the theoretical calculation of Ref. [IJ] with the 
aim of comparing it with the recent JLab data. The spectral function that 
describes the energy and momentum distribution of the nucleons in nuclei is 
obtained by using the Lehmann's representation for the relativistic nucleon 
propagator and nuclear many body theory is used to calculate it for an in- 
teracting Fermi sea in nuclear matter [isl]. A local density approximation is 
then applied to translate these results into finite nuclei. The contributions 
of the pion and rho meson clouds are taken into account in a many body 



field theoretical approach which is based on Refs. [IJ, |16[. The model from 
Ref. [IJ] has been improved in several ways. The old model used the Bjorken 
limit and assumed the Callan-Gross relationship for nuclear structure func- 
tions F2'^{x) and Fi^{x). Due to the fact that JLab data have been taken 
in a region of relatively low (Q^ ~ 3 — 6 GeV^) we have not assumed 
the Bjorken limit. Also, for low and moderate x values Target Mass 
Corrections (TMC) might play an important role. We have incorporated 
them following Ref. [l7|. Another difference with respect to Ref. IJ] is the 
fact that for the ratios we divide by the deuteron structure function, rather 
than the nucleon one. This only implies substantial changes at moderate 
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Figure 1: Feynman diagram for the deep inelastic lepton-nucleon scattering 



and high x values. We have also considered shadowing because it reduces 
the contribution coming from the pion and rho meson clouds 18|, |l9|]. For 



the numerical calculations, next to leading order (NLO) Parton Distribution 
Functions (PDF) for the nucleons have been taken from the parameteriza- 
tion of Martin et al. (MSTW) The NLO evolution of the deep inelastic 



structure functions has been taken from the works of Vermaseren et al. 21 



and van Neerven and Vogt [22]. In the case of pions we have taken the pio- 



nic parton distribution functions given by Gluck et al. 23|, l2J] . For the rho 



mesons, we have applied the same PDFs as for the pions as in Ref. [14 



The structure of the paper is as follows: In Sect. [2] we introduce some 
basic formalism for lepton-nucleon scattering, in Sect. [3] we analyse the dif- 
ferent nuclear effects, in Sect. H] we consider the deuteron case and we end 
by comparing our results with data in Sect. [5l 

2. Deep inelastic lepton-nucleon scattering 

The double differential cross section for the reaction of scattering of a 
charged lepton from an unpolarized nucleon in the one photon exchange 
approximation, 

r (A;) + N{p) l-{k') + X{p'), 1= e, (1) 
depicted in Fig{T]is given, in terms of the Bjorken variables x and y, by 



dxdy 



xy^F,{x,Q^)+{l-y 



xyM\ 



2E 



) 



(2) 
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where 



and V is the energy transferred to the hadronic system. -Fj(x, Q"^) are dimen- 
sionless structure functions. In the Bjorken hmit, i.e. — ?■ oo, v — )■ cxd, x 
finite, the structure functions Fjjx, Q"^) depend only on the variable x and 



satisfy the Callan-Gross relation |25| given by 2xFi{x) = F2{x). Using this, 
the cross section of Eq.(l2]) can be expressed in terms of ^2(0;) and thus the 
ratio of cross sections is equal to the ratio of structure functions F2. Even 
far from the Bjorken limit or when one goes beyond the lowest order (LO), 
where the Callan-Gross relation does not hold, the ratio of cross sections 
still equals the ratio of structure functions F2 if the ratio of longitudinal to 
transverse cross sections R = ai/ does not depend on A. There is a con- 
siderable amount of experimental evidence supporting this fact (e.g. Fig. 6 
of Ref. Therefore, in the following we only consider F2 and compare 

directly F2 ratios with cross section ones. 

The nucleon structure functions are determined in terms of parton distri- 
bution functions for quarks and anti-quarks. In this work, for the nucleons we 
work at NLO and we have used the Parton Distribution Functions (PDF) of 
Martin et al. (MSTW) [i^. At this order, the expression for the ^and Fr, 
structure functions can be expressed as functions of the PDFs by j21, 2^, 2y 



X-^F2,L = ^2,L ® /, (4) 



where C2,l are the coefficient functions for the quarks and gluons [2l|, |22|, |26 



and / represents the quark and gluon distributions [20 



3. Nuclear effects 

We have used the local density approximation (LDA) to incorporate nu- 
clear medium effect^. Inside the nucleus, when the reaction given by Eq.([T]) 
takes place, several nuclear effects like Fermi motion, binding, pion and rho 



^On the other hand, the leading order (LO) pionic parton distribution functions of 
Gluck et al. 23|, |2J] have been used for pions as well as for rho mesons. 
^The nuclear densities have been taken from Ref. 27| . 
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meson cloud contributions must be taken into account. Fermi motion and 
nucleon binding are implemented through the use of a nucleon spectral func- 
tion. The relativistic nucleon propagator in a nuclear medium can be cast 
as [IJ, ^h 



G{p) 



M 



^u,(p) 



du)— : h 



17] 



duj- 



— u + it] 



(5) 



where Sh{oj,p) and Sp{u,p) are the hole and particle spectral functions re- 
spectively. Full details can be found in Ref. [15[. We ensure that the spectral 
function is properly normalized and we get the correct Baryon number for 
the nucleus. Furthermore, we have also calculated the kinetic energy and 
the binding energy per nucleon and have found that the theoretical binding 
energy is very close to the experimentally observed ones for ^Be, ^^C, ^°Ca 
and ^^Fe. 

Our base equation for the nuclear structure function in an isoscalar 
target is: 



d^r 



d^p 



3, dujS,iu,p,pir))^-lMl (6) 



X I 7" + ^ 



— oo 



Pz 



T 



with p° = M + 7'^ = 1 + 4a;'^»VQ^ and x' is Q'^/{2p ■ q). This expres- 
sion is equivalent to that of Ref. 18| after trivial algebraic transformations 
and taking into account the different normalization of the spectral function 
Vo{e,p) used in [18] such that 



A Voie, p) 4 ■ 27r y c/V Shitu, p, p(r)) . 



(7) 



In an earlier study, the behaviour of different nucleon spectral func- 

In particular, the spectral functions given by 



tions has been analysed [28 



Fernandez de Cordoba and Oset [15| , Kulagin and Petti [18| , and Ankowski et 
al. [29| were used and compared. It was found that the results do not change 
appreciably. Finally, we should comment that the present formalism has also 
been used to study the nuclear effects in the F3 structure function [30 . 
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3.1. vr and p mesons contribution to the nuclear structure function 

The pion and rho meson cloud contributions to the F2 structure function 
have been iniplemented following the many body field theoretical approach 



of Refs. 2M, 116|. The pion structure function F2a,tt{x) is written as 



F^^ix) = -Q [ d?r [ 4^ ^(/) 6ImD{p) — 2M F^^ix^) e{x^-x) e{l-x^) 

(8) 

where D{p) the pion propagator in the medium given in terms of the pion 
self energy 11^: 

Dip) = _ _ ^2 _ n,(p0^p)]-i ^ 

where 

P/mlF\p)p'U* 
" l-p/mlV[U* ■ ^ ^ 

Here, F{p) = (A^ - ml)/{K^ + P ^) is the vrA^A^ form factor and A=l GeV, 
/ = 1.01, V[ is the longitudinal part of the spin-isospin interaction and 11* 
is the irreducible pion self energy that contains the contribution of particle - 
hole and delta - hole excitations. In Eq.(IH]), SImD{p) is given by 

6ImD{p) ^ ImD{p) - p ^IH^M (n) 



and 



x^ M ^ ^ 



Assuming SU(3) symmetry and following the same notation as in Ref. j23|. 
the pion structure function at LO can be written in terms of pionic PDFs as 

/5 12 \ 

F2niXn) = X^i- V^{x^) + — q^{x^) j (13) 

where Vn{xn) is the valence distribution and qnix^) is the light SU(3)-symmetric 
sea distribution. 

Similarly, the contribution of the p-meson cloud to the structure function 
is written as 14 
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/f 
d'r J -^^e{p'')5ImDp{p)^2MF2p{xp)e{xp-x)e{l-Xp) 

(14) 

where Dp{p) is the p-meson propagator and F2p{xp) is the p- meson structure 
function, which we have taken equal to the pion structure function F2-K using 
the valence and sea pionic PDFs from reference j23|. Kp in pNN form factor 
F{p) = (A^ — m?p)/{JS?p + P^) has also been taken as 1 GeV. 

Further details concerning the pion and p-meson propagator can be found 
in Ref. [l^. This model for the pion and p selfenergies has been abundantly 
used in the intermediate energy region and provides a quite solid description 
of a wide range of phenomenology in pion, electron and photon induced 



reactions in nuclei, see e.g. Refs. [3l|, |32|, |33|, [SJ, l35[ and references in [14 



In particular, a careful study of the in medium pion propagator used here 



was carried out in Ref. [16|. There, several tests concerning the fulfillment 
of sum rules, and the preservation of the analytical properties of the meson 
propagator and the consistency of the results with similar calculations were 
considered. 

In addition, the balance of light-cone momentum between bound nucleons 
and pions can be studied by means of a momentum sum rule as done in 
Ref. [l8]. The pion < ?/ >^ and nucleon < y >n fractions of the light cone 
momentum are related by 

<y>^ + <y>N=j^, (15) 



where Ma is the nucleus mass. See section 5.3 of Ref. [18| for details. The 
sum rule should be valid for a nuclear model where the Hamiltonian would 
contain only pions and nucleons. In fact, our model for the nucleon spectral 
function is based on a phenomenological approach that also contains many 
other pieces in the nucleon-nucleon interaction and thus the sum rule is not 
directly applicable. Nonetheless, it can provide further constraints on the size 
of the mesonic contribution and it will be discussed in the results section. 

The mesonic cloud contribution is expected to be negligible for deuteron 
as it depends, roughly speaking, quadratically on the baryon density which 
is quite small for this cas^. Therefore, these contributions have not been 
included in the evaluation of the deuteron structure function. 



A direct application of our model to deuteron, produces a mesonic contribution that 
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3.2. Target mass corrections 

Target mass corrections have been incorporated by means of the approx- 
imate formula 17 



?TMC , 



X 



^2 y 



1 + ^(1-0^' 

7 



(16) 



where /i = 7 



1 + ^^^y^ and ^ is the Natchmann variable defined as 



2x 

1 + 7 



(17) 



3.3. Coherent nuclear effects 

Furthermore, we have taken into account the shadowing effect following 
the works of Kulagin and Petti 18[ . We are interested in the relative effect 
in F2 that can be written as 



5Ro 



5F^ 1 + 



1 + R 



5Ri 



(18) 



where R{x, Q^) is calculated for the free nucleon. For 5Rt we use expression 
63 in Ref. |l8|. The most general expression for i?(x, Q^), taking into account 
the target mass is: 



R{x,Q^) 



Fl 7'^2 - 2a;Fi 72^2 



Ft 



2xFi 



2xFi 



- 1. 



(19) 



4. Derivative expansion of F2 

The difference between dividing by the deuteron F,p or by the free isoscalar 
structure function to calculate the ratios of structure functions is of only 
a few percent in the x < 0.7 region. However, the quality of data requires 
a proper description of F^. On the other hand, our local density approach 
is not appropriate for such a light nucleus (or even for ^He). Therefore, we 



is always lower than a 0.6 percent of the nucleonic contribution for the analysed x range. 
Thus, its inclusion would have a very minor effect in the ratios. Nonetheless, we should 
mention that our formalism, which starts from selfenergies calculated in nuclear matter, 
is not expected to be very reliable for the calculation of the mesonic effects in deuteron. 
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need another method to calculate F2 in these cases. In terms of the deuteron 
wave function, can be written as 



(20) 

Alternatively, a particularly appealing approach because of its simplicity, 
is the use of derivative expansions that provide the structure function per 
nucleon of a nucleus in terms of the free nucleon one, its derivatives and a 
few expected values of nuclear observables 28, 37, 38, 39, 4^. We can write 



where < T > is the mean nucleon kinetic energy taken as 11.07 MeV, < E > 
is the nucleon removal energy taken as 2.226 MeV, < Tr >~ (p^) /2M with 
(p^) = 0.533 /m~^ the average of the square of the nucleon momentum. 
To include TMC, one must substitute in Eq. ( !2T|) the free nucleon structure 
function and its derivatives by the approximate one given in Eq. ( |T6l) . 

The derivative expansions have some intrinsic limitations and it has been 
shown that they fail to converge to the results obtained by folding with the 
nuclear spectral functions for x > 0.5, for the case of medium and heavy 
nuclei. A detailed study can be found in Ref. ^8|. The convergence is 
expected to be much better for a loosely bound nucleus such as the deuteron. 



Indeed, obtained using the Paris wave function 411] and the results of the 
derivative expansion differ by less than 0.6 percent up to x = 0.6 as shown 
in Fig. H 



In Fig. [21 we also include for comparison the same ratio from Ref. 42 
which uses a different set of PDF's. We have observed that the main differ- 
ence with that calculation comes from the inclusion of a parametrization of 
the off-shell effects (see dashed-dotted line) absent in our model. 



5. Results and discussion 

Our aim in this paper is to confront the model with the recent JLab 
results of Ref. [l| that correspond to ratios of nuclei with deuteron and more 
precisely with the slope of the x dependence that is more insensitive to the 
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Figure 2: F^/{F^ + F^) as a function of x at Q^^IO GeV^ Solid line: Derivative 
expansion. Dashed line: Eq.[20l Dashed dotted: Eg. [20l including ofF-shell effects following 
the prescription of Ref. [18j. Dotted line: Ref. (4^ . 
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Figure 3: Electromagnetic structure function in deuteron F^ at different x values. The 
theoretical curves are obtained by using Eq. (|20p . Experimental data are taken from 
Ref. [13. 
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Figure 4: Electromagnetic structure function in deuteron at different x values. Solid 
lines are obtained by using Eq. (pn|) and the dotted lines are the results for the free nucleon 
case. Experimental data are taken from Ref. 44 1. 



normalization uncertainties. Nonetheless, we will also show some results for 
the deuteron FP(x, Q^) structure function as well as the ratio R(x,Q^) = -|^ 

in intermediate mass nuclei like ^°Ca and ^^Fe. 

In Fig. [31 we compare the theoretical calculation obtained using Eq. ( 120|) 
with the experimental results of Benvenuti et al. jisj. Overall, we find a 
good agreement in the x region relevant for our study. Although the data 
correspond to large values, this gives us confidence in the quality of this 
approach for the evaluation of the ratios with respect to other nuclei. In 
Fig. m we compare our results with data obtained with a muon beam on a 
deuterium target 4^. We also show the results for the free nucleon case. 
The nuclear corrections are very small for the range of x values analysed in 
the experiment. 

One of the most interesting results of the recent JLab data is that for both 
Beryllium and Carbon the cross section ratios show a similar slope even when 
they have a quite different average nuclear density. This conflicts with some 
simple fits that describe well the slope for medium and heavy nuclei as a 
function of the average nuclear density or with simple A dependences jsj. 
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Figure 5: Ratio R(x,Q )= . Full model (solid line), without shadowing (dotted line) 

and without pion, rho and shadowing contributions (dashed line). For each value of a:, 
has been calculated using an electron beam of 5.767GeV and scattering angle of 40° 
corresponding to JLab kinematics. Data are cross section ratios from Ref. 

On the other hand, the slope of the ratio in the region 0.3 < x < 0.6 is 
particularly well suited to analysis because from the experimental point of 
view it is quite unaffected by normalization uncertainties. Also theoretically 
it is relatively simple because shadowing, or Fermi motion are of a little 
importance over this region of x. 

In Fig. [5l we show the results for Beryllium. The dashed line has been 
calculated using Eq. ([6]) with TMC and the solid line corresponds to the full 
model, including the meson cloud contributions, shadowing and TMC. We 
show explicitly the effect of shadowing. It reduces the structure function 
ratio by around 1% at a; ~ 0.3 and even less for higher x. We have found 
that TMC has a really minor effect in the ratio for these x values (less than 
1% at a; ~ 0.6 and even smaller for lower x values). Therefore, the difference 
between the base curve and the full one comes basically from the vr and 
p contributions that play an important role. The size of the rho meson 
correction is about half that of the pion. We find that the full model agrees 
quite well with data both in slope and the size of the ratio. 

A good agreement with data is also obtained for Carbon as shown in 
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Figure 6: Ratio R(x,Q^)= ^^^^ ■ Full model: solid line with A, Ap ~ IGeV; the band 

corresponds to ±20% variation on A and Ap. Full model without pion, rho and shadowing: 
dashed line. for calculation and JLab data [1] as described in previous figure. SLAC 
data correspond to = 5 GeV^ . 



Fig. [6l The slope and size of the nuclear effects are similar to the Beryllium 
case. This could look surprising given the quite different average density as 
discussed in [l|. This points out to the fact that "average density" could not 
be the appropriate parameter for the description of the EMC effect in light 
nuclei. For example, this has been discussed in Sect. IV of Ref. jsol]. Again, a 
determining factor in the agreement is the mesonic cloud contribution. Given 
this, some words of caution are needed here. First, the parton distribution 
functions are poorly known for the mesons and possible off-shell effects have 
not been included in the calculation. Second, the results depend on the 
meson selfenergies in the medium that also contain some uncertainties such 
as the specific form of the spin-isospin interaction, specially for the p meson. 
A full analysis of these uncertainties is out of the scope of this paper. To give 
an idea of their size, we have shown in this figure the results for the ratio 
using the full model with A, Ap=lGeV and A, Ap=1.2GeV and O.SGeV. We 
find that a 20% variation in the A's, results in a 2-3% change in the ratio. 

We have also tested the momentum sum rule discussed in section 13.11 
The mesons carry a light-cone momentum fraction of 3 percent for A = 1 
GeV. The 2 percent prescribed by the sum rule can be obtained for a cut-off 
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A = 0.8 GeV. This suggests that lower cut-off values should be preferred 
but one must be careful before reaching such conclusion. For instance, the 
nucleon momentum fraction is very sensitive to parameters like the expected 
value of the nucleon kinetic energy that are not very well known and has some 
uncertainty. We have used the values obtained with our nucleon spectral 
function. The same results are obtained for ^Be. Heavier nuclei, such as iron 
and calcium, have a mesonic momentum fraction of 5 percent and fulfill the 
sum rule for A = 0.75 GeV. 

In Fig. [6l the systematic difference in size between JLab and SLAG data 
is consistent with the normalization uncertainties quoted in Refs. BQ. It 
may be noted, however, that the slope is very similar for both experiments 
and in good agreement with our results. These normalization differences 
have been recently discussed in Ref. [42|. 

For both nuclei, our results slightly overestimate data by around 2% at 
X around 0.6 and more above that. However, that region is much affected 
by possible off-shell effects [42], not included in our approach, and by high 
momentum components of the nucleons wave function. Therefore, we cannot 
make any strong statement about this discrepancy apart from the fact that 
we are reaching one of the limits of validity of our model. 

We have also checked that the use of next to next to leading order PDFs, 
that considerably lengthens and complicates the calculation, does not appre- 
ciably change the results, at the level of precision of the current data and the 
size of other theoretical uncertainties. 

There are JLab results even for lighter nuclei like ^He and ^He. Our local 
density model is certainly not adequate for these cases that would require a 
more microscopical approach for the calculation of a proper nucleon spectral 
function and of the meson cloud contribution. Also, good data for larger 
X values are available. They are particularly sensitive to TMG and to high 
momentum components of the nucleon spectral functions. In order to analyse 
these data, further work would be required to extend the validity of the 
theoretical approach describing the nucleon spectral function. 

As a further test, we have also studied the results for the ratio R(x,Q^)=^ 
for intermediate mass nuclei like calcium and iron. The results are shown 
in Figs. [Tland [HI In both cases, we have compared with SLAG results 
from Ref. |5[, with averaged Q^. No significant dependence was found 
in Ref. jsf for their kinematic range. The theoretical curves have been cal- 
culated for = 5GeV^ and we observe little sensitivity to that value. In 
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Figure 7: Ratio R(x,Q^)= ^p . Full model: solid line. Full model without pion, rho 
and shadowing: dashed line. Dotted curve is the full result scaled by a factor 0.97. 
Calculations have been done for = 5 GeV^. The experimental points are taken from 
Ref. (averaged Q^). 



the case of calcium, our results overestimate the data by around 3%. This is 
larger than the normalization uncertainties quoted in Ref. jsj. Nonetheless, 
scaling our theoretical curve, we observe a good agreement with the slope of 
the structure function. 

The situation is much the same for heavier nuclei, such as iron. The slope 
is well reproduced and calculation overestimates again data by around a 3%. 
Similar results are obtained for silver and gold. This overestimation seems to 
be consistent with the results of recent global fits to the nuclear parton distri- 



bution functions (see e.g. Fig. 4 of Ref. [13|). Their results might point out 
to some normalization uncertainty in the SLAC results such that the medium 
and heavy nuclei ratios are too small. However, the recent more microscopi- 
cal analysis of Ref. 42| that also finds normalization inconsistencies between 
light and heavy nuclei favours the interpretation that the recent JLab data 
should be rescaled by a global factor of around 0.98 and that the SLAC data 
are correct. We should also mention that the two discussed experiments had 
significantly different lepton energies and the simple A dependence assumed 
for the ai/ax ratio could have to be revised. 

Certainly, these normalization issues should be settled with new and bet- 
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ter experiments. From the theoretical point of view, it seems that microscop- 
ical models are hardly able to reproduce at the same time the high statistic 
data from light nuclei at JLab and medium and heavy nuclei from other 
collaborations. In any case, this does not affect the main point discussed 
in this paper, namely the slope produced by the nuclear effects that is well 
reproduced in our model. 

In summary, the electromagnetic nuclear structure function F2 has been 
studied including nucleonic and mesonic degrees of freedom for a x region 
where shadowing, antishadowing and Fermi motion are not too important. 
We have started from up to date nucleonic PDF. Nuclear effects like Fermi 
motion and binding have been incorporated by means of the use of a spectral 
function obtained for nuclear matter and implemented in nuclei using the 
local density approximation. A similar approach has been used for the inclu- 
sion of the contribution of mesonic clouds. Also shadowing and TMC have 
been considered. The deuteron structure function has been calculated using 
a derivative expansion and with the Paris wave function. The results success- 
fully reproduce recent very precise JLab results for light nuclei at relatively 
low values. Also the slope of previous experiments for heavier nuclei is 
well reproduced although we fail to agree with them on the absolute size by 
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up to a 3%, larger than the quoted experimental uncertainty. We have found 
that the mesonic cloud (basically pion) gives an important contribution to 
the cross section ratios but it still has considerable uncertainties. Even small 
changes of the pion nuclear selfenergy can produce appreciable changes in 
the cross section ratios. 

The success of this local density model for light nuclei is in contrast with 
the failure of simple models/parametrizations that fit well for the nuclear 
effects for medium and heavy nuclei as a function of average density or the 
mass number A The use of an approach that incorporates in an adequate 
manner the nucleon and meson properties in the nuclei is clearly mandatory 
for the analysis of the EMC effect in these cases. 
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